Update on multiparticle effects in 
Bose-Einstein correlations 



Kacper Zalewski* 
M. Smoluchowski Institute of Physics, Jagellonian University 

and 

Institute of Nuclear Physics, Krakow, Poland 
February 2, 2008 



Abstract 

Multiparticle effects in Bose-Einstein correlations are reviewed. It 
is shown that for a broad class of models they can be ignored in the 
low density limit, but often are significant (typically at a 10% level) 
for realistic final states. 

PACS 25.75. Gz, 13.65.+i 
Bose-Einstein correlations. 



1 Statement of the problem 

Most models of the Bose-Einstein correlations start with an input single par- 
ticle density matrix p/(p;p') or some equivalent information. The subscript 
/ reminds that this is an input and not necessarily the actual single parti- 
cle density matrix of the system being considered. This input depends on 
the quantities of interest e.g. on the radii of the interaction region. Let 
us concentrate on a system of n identical bosons. If these particles were 
distinguishable, one could build a reasonable n-particle density matrix as a 
product of the single particle density matrices pf. 
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Pn(Pl,---,Pn) = I[Pl{PjiP'j)- C 1 ) 

i=i 

For identical bosons, however, this product must at least be symmetrized pQ , 
[2] and one finds 

n 

Ai(pi,---,Pn) = SIIMPi;ppi)» ( 2 ) 

P j=l 

where the summation is over all the permutations of the momenta p' 1; . . . , p' n . 
The tilde reminds that in general Trp ^ 1 and thus, strictly speaking, p n is 
not a density matrix. This matrix is generally used to calculate the /c-particle 
momentum distributions Qk(P\, ■ ■ ■ , Pjt) putting n = k. Thus 



Q 1 (p) = C lPl (p;p), (3) 

^2(pi,p 2 ) = C 2 [p/(p 1 ;pi)p/(p 2 ;p2) 
+ p/(pi;p 2 )p(p 2 ;pi)] , ( 4 ) 

etc., where Cj are normalization constants. This is an approximation. 
should have contributions from all n > k. A way of taking that into account 
has been proposed by Pratt Suppose that without symmetrization the 
probability of producing exactly n identical bosons is Po(n) - this is another 
input. Then it is natural to assume that with symmetrization this probability 
goes over into 

Pin) = C P {n)Trp n , (5) 
where by definition Trp = 1. Formulae (jBJ) etc. get replaced by 

^i(p) = p/(p;p) + 

oo „ 

X)-Po(n) / dp 2 ,...,dp„/o ri (p 1 ,...,p n ;p lJ ...,p n ), (6) 

n=2 

etc., where the normalization constants have been omitted. Depending on 
conventions, dp may mean either d 3 p or The difference between formulae 
(J3j) and analogous and the corresponding formulae (0) and analogous is due to 
the terms with integrations, or in other words to the multiparticle effects. The 
problem is (cf. e.g. [1], how important in practice are these multiparticle 
effects? 
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2 Low density limit 



Since every permutation can be decomposed into cycles, all the integrals 
occurring in Pratt's model can be expressed in terms of two types of basic 
integrals. These are the chain integrals: G\(p; p') = pj(p; p') and for m > 1 



Gm(p,p') = 

/ p/(p; Pi)dpiP/(Pi; P2) • • • dpm-i/oj(p m -i; p'); ( 7 ) 

and the cycle integrals 

C m = J dp G m (p;p). (8) 

According to this definition C\ = 1. The evaluation of the basic integrals 
can be greatly simplified [Zj when the eigenfunctions and eigenvalues of 
the density matrix pj, 

J Pi(p;p')dp'ipk(p') = Afe^(p), (9) 

are known. Then 



G m (p;p') = p?(p;p') = 

E^(p)W(p), (10) 

k 

C m = Trp? = J2W- (11) 

k 

All the eigenvalues of a density matrix are real, nonnegative and not exceed- 
ing one. Let us label them so that 



< ... < Ai < A Q < 1. (12) 

The most populated state has the occupation probability (per particle) A . 
Therefore, the low density limit corresponds to A —>■ 0. In this limit all the 
basic integrals except C\ and G\ tend to zero and, consequently, all the inte- 
grations drop out. In the low density limit formula (JHJ) and its more-particle 
analogues become exact - the multiparticle effects can be ignored. The limit 
of Ao is taken at constant Po{n). This corresponds to a constant average 
number of particles before symmetrization and, therefore, using the classical 
terminology, the density tends to zero, because the amount of available phase 
space grows to infinity. 
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Let us consider an example. Pratt [3] chose 



P/(p; P') = ; 3 e-^*—* R (pp } . (13) 



V27TA2 

For this choice 0, 0, |H1 



A ° = [tTIKr) ■ (14 > 

The parameter A 2 (R 2 ) is equal to the average of the square of a component 
of the momentum vector (position vector). Thus, in the limit Ao — > 0, which 
corresponds to AR — ► oo, the phase space available tends indeed to infinity. 



3 Finite density 

When A does not tend to zero the results become much more complicated 
. The two-particle correlation function, which according to (J3J) is 

C 2(Pi,P 2 ) = H 7 \ 7 v ( 15 ) 

p/(pi;pi)p/(p 2 ;p 2 ) 

becomes 

c 2 ( Pl , P2 ) = cr(Pi, p 2 ) [i + iMPi, p 2 )] > (is) 

where 



rres , s _ E,,^(^)^(Pl.Pl)^(P 2 P 2 ) n? x 

<- 2 lPl,P 2 J - i «,(t)«,(,-)Gf 1 (p li p 1 )Oi(P 9 .P a )' 1 J 

«2lPi,P 2 J - X) iii «(i+i)Oi(P 1 ,P 1 )Oj(P 2 ,P a )' 

^M=E~ m P(n)^^. (19) 

Function i? 2 is qualitatively similar to the second term on the right-hand 
side of formula (j!5)l . It is equal one for p x = p 2 and tends to zero for 
( Pl — p 2 ) 2 — > oo. It can be used to extract information about the geometry 
of the interaction region. Function C 2 es has been called residual correlation 
0> [O]- We will see further that it is identically equal one when the input 
multiplicity distribution Po(n) is Poissonian, but in general it complicates 
significantly the analysis. In particular, it can make the value of C 2 (p 1; p 2 ) 
at large (p 1 — p 2 ) 2 bigger than at p 1 = p 2 jS]. 
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Let us consider as an example the case when the input multiplicity dis- 
tribution is Poissonian with average multiplicity v: 

P (n) = — s e- v . (20) 
n! 

Then w(m) = z/ m , which implies w(i + j) = w(i)w(j), and consequently 
C2 £S (Pi)P2) = 1- Moreover, introducing the notation p] 

oo 

L(p;p') = E^^(p;p'), (21) 

3=0 

one finds [7j 

g a (Pi,p a )=i+ 7? L(P \ ;P ; )|2 , (22) 

^(pi;pi)^(p 2 ;p 2 ) 

which is formula with the replacement pi — > L. Similar relations hold for 
more-particle distributions and correlations. This analogy, however, can be 
misleading. Interpreting L(p 1 ;p 2 ) as if it were /9j(p 1 ;p 2 ) one overestimates 
the average particle multiplicity (n) and underestimates the averages (p 2 ) 
and (x 2 ). Thus, the phase space particle density is overestimated. 



4 How important is all that in practice ? 

The multiparticle effects are negligible when the average population of each 
state, or using classical terminology the particle density in phase space, is 
sufficiently low. They are very important when this population, or density, 
is very high. The obvious question is: what about the actual data ? We will 
present two kinds of estimates suggesting that the multiparticle corrections 
are likely to be of the order of 10 %. 

4.1 Monte Carlo approach 

Most Monte Carlo codes first generate events without Bose-Einstein corre- 
lations. These correlations are introduced later by modifying the weights of 
the events or by shifting the momenta of the generated particles. The for- 
mer approach has a problem: the uncorrected Monte Carlo probabilities give 
multiplicity distributions in agreement with experiment; the weights modify 
the multiplicity distribution and thus tend to spoil its agreement with exper- 
iment. A way out is (cf. e.g. jTU], [TT] ) to introduce two new parameters: c 
and V satisfying approximately the relations 



5 



P (n) = cV- n P w (n). (23) 

Here Po(n) is the uncorrected multiplicity distribution which, by construc- 
tion, agrees with experiment. P w {n) is the distribution modified by the 
weight factors. A priori there is no reason why two parameters should be 
enough to get reasonable agreement for the full multiplicity distributions, but 
in practice it works. The deviation of the parameter V from unity is a mea- 
sure of the importance of symmetrization effects. Very roughly, it gives the 
change of the average particle multiplicity as result of the symmetrization. 
It was found that V = 1.1 for pp scattering at centre of mass energy ^/s = 30 
GeV HU and V = 1.05 for e + e~ annihilations at = 172 GeV [TU|. Thus 
the effects are small, but not negligible when a precision of some 10 % is 
aimed at. These results are model dependent. It is possible to introduce 
weights which oscillate around one so that V ~ 1 J2j, jE] Symmetrization 
does not affect the multiplicity distribution at all in the method of momen- 
tum shifts jHj, ^21- The shifts affect neither the weight of an event nor its 
multiplicity. In this model, however, the authors also find [Tj)| that multipar- 
ticle effects are important, though not for the multiplicity distribution which 
is protected by construction. One important effect of symmetrization, com- 
mon to many models, is the enhancement of the single particle momentum 
spectra at low momenta. 

4.2 Thermodynamics 

Many models, not only the models called thermodynamic, introduce the 
temperature T and the chemical potential \x of the particles. In such models 
the momentum spectrum of the particles is given before symmetrization by 
the Maxwell-Boltzmann formula 

Omb(p) = eTe~^ (24) 
and after symmetrization by the Bose-Einstein formula 

<Wp) = • (25) 

e t — 1 

The proportionality constant in the Maxwell-Boltzmann formula has been 
chosen so that at large values of the particle energy E p the formula coincides 
with the Bose-Einstein one. The effect of symmetrization is the "—1" in 
the denominator of the Bose-Einstein formula. The first remark is that this 
correction becomes more important with respect to the exponential when 
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temperature increases. Incidentally, for this reason some people consider 
Einstein's condensation a high temperature effect. 

Temperatures in high energy multiple particle production processes are 
strongly model dependent. Heavy ion collisions models where all transverse 
momentum is due to thermal motion have T ~ (pt) ~ 350 MeV. Models 
where much of the transverse momenta is due to collective radial expansion 
have lower (local) temperatures of about 100 MeV. The chemical potential 
cannot exceed the lowest possible energy, in this case m n = 140 MeV. Fits 
assuming temperatures of about 100 MeV [TH], [Ej give values of \x ranging 
from 30 MeV to 90 MeV. We will give two estimates of the importance of 
multiparticle effects based on these numbers. 

The symmetrization effects are the strongest at p = 0. The ratio of the 
Bose-Einstein to the Maxwel-Boltznann distribution at p = is 

n = K^r - (26) 

1 — e t 

For T = 100 MeV and /i ranging from 30 MeV to 90 MeV r\ ranges from 
1.5 to 2.5. Thus the correction is very significant. A more representative 
ratio is the ratio of the integrated multiplicities corresponding to the two 
distributions: 

Jp 2 Q BE {p)dp 

r2 = T^ri TXT- ( 27 ) 

jp 2 VL MB {p)dp 

For the same choice of parameters r 2 ranges from 1.09 to 1.20, which is 
consistent with the Monte Carlo results. 



5 Conclusions 

The importance of the multiparticle effects depends both on the model be- 
ing used and on the quantities of interest. In Monte Carlo models, fine 
tuned to reproduce without symmetrization the multiplicity distributions, 
symmetrization is forbidden to affect the multiplicity distributions. In mod- 
els related to thermodynamics the low momentum particles are affected much 
more than the high momentum particles. It is not unusual to get from multi- 
particle effects corrections of the order of 10 % or more. Therefore, it seems 
that in precision work the multiparticle effects should not be discarded with- 
out previous analysis. 
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